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We report on the first realization of a single bosonic Josephson junction, implemented by two 
weakly linked Bose-Einstein condensates in a double-well potential. In order to fully investigate 
the nonlinear tunneling dynamics we measure the density distribution in situ and deduce the evo- 
lution of the relative phase between the two condensates from interference fringes. Our results 
verify the predicted nonlinear generalization of tunneling oscillations in superconducting and super- 
fluid Josephson junctions. Additionally we confirm a novel nonlinear effect known as macroscopic 
quantum self-trapping, which leads to the inhibition of large amplitude tunneling oscillations. 
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Tunneling through a barrier is a paradigm of quan- 
tum mechanics and usually takes place on a nanoscopic 
scale. A well known phenomenon based on tunneling is 
the Josephson effect [l( between two macroscopic phase 
coherent wave functions. This effect has been observed 
in different systems such as two superconductors sepa- 
rated by a thin insulator 0] and two reservoirs of super- 
fluid Helium connected by nanoscopic apertures (3, La]- 
In this letter we report on the first successful implemen- 
tation of a bosonic Josephson junction consisting of two 
weakly coupled Bose-Einstein condensates in a macro- 
scopic double-well potential. 

In contrast to all hitherto realized Josephson junctions 
in superconductors and superfluids, in this new system 
the interaction between the tunneling particles plays a 
crucial role. This nonlinearity gives rise to new dynami- 
cal regimes. Anharmonic Josephson oscillations are pre- 
dicted d, IE ! if the initial population imbalance of the 
two wells is below a critical value. The dynamics changes 
drastically for initial population differences above the 
threshold of macroscopic quantum self-trapping || 0, 
where large amplitude Josephson oscillations arc inhib- 
ited. The two different dynamical regimes have been 
experimentally inve stig ated in the context of Josephson 
junction arrays [ill Il2l fl3 | . However, the small period- 
icity of the optical lattice does not allow to resolve indi- 
vidual wells and thus the dynamics between neighboring 
sites. Our experimental implementation of a single weak 
link makes it possible for the first time to directly ob- 
serve the density distribution of the tunneling particles 
in situ. Furthermore we measure the evolution of the 
relative quantum mechanical phase between both con- 
densates by means of interference 0] . 

The experimentally observed time evolution of the 
atomic density distribution in a symmetric bosonic 
Josephson junction is shown in Fig.^for two different ini- 
tial population imbalances (depicted in the top graphs). 
In Fig. IHa) the initial population difference between the 
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FIG. 1: Observation of the tunneling dynamics of two weakly 
linked Bose-Einstein condensates in a symmetric double- well 
potential as indicated in the schematics. The time evolution 
of the population of the left and right potential well is directly 
visible in the absorption images (19.4 /jm x 10.2/nn). The 
distance between the two wavepackets is increased to 6.7/im 
for imaging (see text), (a) Josephson oscillations are observed 
when the initial population difference is chosen to be below 
the critical value zc- (b) In the case of an initial population 
difference greater than the critical value the population in 
the potential minima is nearly stationary. This phenomenon 
is known as macroscopic quantum self-trapping. 

two wells is chosen to be well below the self-trapping 
threshold. Clearly nonlinear Josephson oscillations are 
observed i.e. the atoms tunnel right and left over time. 
The period of the observed oscillation is 40(2)ms which 
is much shorter than the tunneling period of approxi- 
mately 500ms expected for non-interacting atoms in the 
realized potential. This reveals the important role of 
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the atom-atom interaction in Josephson junction experi- 
ments with Bose-Einstein condensates. A different man- 
ifestation of the nonlinearity is shown in Fig. ^b) ex- 
hibiting macroscopic quantum self-trapping, which im- 
plies that the population imbalance does not change over 
time within the experimental errors. The only difference 
to the experiment shown in Fig. ^a) is that the initial 
population imbalance is above the critical value. 

The experimental setup and procedure to create the 
87 Rb Bose-Einstein condensates is similar to that used 
in our previous work 0|. A sufficiently precooled ther- 
mal cloud is loaded into an optical dipole trap consisting 
of two crossed, focussed laser beams and is subsequently 
evaporatively cooled by lowering the light intensities. 
We produce pure condensates consisting of 1150 ± 150 
atoms and final trap frequencies of u> x = 2ir x 78(l)Hz, 
uiy = 2ir x 66(l)Hz and lo z = 2ir x 90(l)Hz, with gravity 
acting in the y-direction. Subsequently we adiabatically 
ramp up a periodic one-dimensional light shift potential 
in x-direction to a depth of 2tt x 412(20)Hz with peri- 
odicity 5.2(2)/im realized by a pair of laser beams at 
a wavelength of 811nm crossing at a relative angle of 
9° . The superposition of this periodic potential with the 
strong harmonic confinement creates an effective double- 
well potential in x-direction with a barrier height of 
2ir x 263(20)Hz, which splits the initial condensate into 
two parts separated by 4.4(2)/^m realizing a single weak 
link (see schematics in Fig. 

The initial population difference between the left and 
right component is obtained by loading the Bose-Einstein 
condensate into an asymmetric double-well potential, 
which is created by a displacement of the harmonic con- 
finement with respect to the periodic potential. The 
asymmetry can be adjusted by shifting the focussed laser 
beam which realizes the harmonic confinement in x- 
direction. This is done by means of a piezo actuated 
mirror mount. A relative shift of only 350nm leads to a 
relative population difference corresponding to the self- 
trapping threshold. This demands high passive stability 
of the mechanical setup and makes it necessary to ac- 
tively stabilize the phase of the periodic potential. With 
our setup we can adjust any initial population imbalance 
with a standard deviation of Az — 0.06. The Joseph- 
son dynamics is initiated at t = by non-adiabatically 
(with respect to the tunneling dynamics) changing the 
potential to a symmetric double-well (see schematics in 
Fig. |2J). After a variable evolution time the potential 
barrier is suddenly ramped up and the harmonic poten- 
tial in x-direction is switched off. This results in dipole 
oscillations of the atomic clouds around two neighbor- 
ing minima of the periodic potential. Thus by imag- 
ing at the time of maximum separation (1.5ms) we can 
observe clearly distinct wave packets with a distance of 
6.7(5)/im. The atomic density is deduced from absorp- 
tion images with a spatial resolution of 2.2(2)/^m. In 
previously reported realizations of Bose-Einstein conden- 



sates in double-well potentials |16l Il7j the time scale of 
tunneling dynamics was in the range of thousands of sec- 
onds. In contrast, our small interwell distance combined 
with a negligible thermal atomic fraction allows the re- 
alization of tunneling times on the order of 50ms, which 
makes the direct observation of the nonlinear dynamics 
in a single bosonic Josephson junction possible for the 
first time. 

The physics of Josephson junctions is based on the 
presence of two weakly coupled macroscopic wave func- 
tions separated by a thin potential barrier. Insight into 
the dynamics of the system can be gained by employing 
a two mode approximation which characterizes the wave 
function by only two parameters, the fractional relative 
population z — (Ni — N r )/(Ni + N r ) and the quantum 
phase difference <\> = <fi r — <f>i between the left (I) and 
right (r) component. In this framework the resulting 
quantum dynamics in a symmetric double-well potential 
is described by two coupled differential equations 

i = — \J \ — z 2 sin0 (1) 
6 = Az-\ — , cos 6 

VT^l 2 

where A is proportional to the ratio of the on-site in- 
teraction energy and the coupling matrix element given 
in ||. These equations represent the nonlinear general- 
ization of the sinusoidal Josephson oscillations occurring 
in superconducting junctions. An intuitive understand- 
ing of the rich dynamics of this system can be gained by 
considering a descriptive mechanical analog. The equa- 
tions given above describe a classical non-rigid pendulum 
of tilt angle </>, angular momentum z, and a length pro- 
portional to y/1 — z 2 . In the following discussion we will 
only consider the case of vanishing initial phase difference 
(f)(0) = 0. If the initial population imbalance is below the 
critical value [Io| |z(0)| < zc (from our experimental 
results we deduce zc ~ 0.5 corresponding to A m 15), 
equ. n describes oscillations in z and <j). In the limit of 
|z(0)| <C zc this reduces to a harmonically oscillating 
mathematical pendulum. In the context of Josephson 
junctions this behavior is known as plasma oscillations. 
A different dynamical phenomenon arises if the initial 
population imbalance is above the critical value. This 
implies that the difference between the two on-site inter- 
action en ergi es becomes larger than the tunneling energy 
splitting |l8j . In this case the relative phase rapidly in- 
creases in time leading to a rapidly alternating tunneling 
current according to equ. ^ This results in a population 
imbalance which performs small oscillations around the 
initial value (self-trapping, running phase modes 0). 

In this case the population difference is self-locked 
to the initial value and the relative phase is increasing 
monotonically (running phase modes |lCJ). In the me- 
chanical analogue this critical imbalance corresponds to 
an initial angular momentum sufficiently large that the 
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FIG. 2: Detailed analysis of the time dependence of the two dynamical variables z and <f> describing the system. The top 
graphs depict the experimental preparation scheme implemented to realize different initial atomic distributions. The dynamics 
is initiated at t — by switching non-adiabatically to the symmetric double-well potential. Graph (a) shows the familiar 
oscillating behavior of both the population imbalance and the relative phase in the Josephson regime. The solid lines represent 
the results obtained by numerically integrating the non-polynomial Schrodinger equation, and are in excellent agreement with 
our experimental findings. The shaded region shows the theoretically expected scattering of the data due to the uncertainties 
of the initial parameters and broadens for large evolution times due to different oscillation frequencies for different initial 
population imbalances. The insets depict representative atomic interference patterns obtained by integrating the absorption 
images along the y- and z-direction after the indicated evolution times. In graph (b) the totally different dynamics in the 
regime of macroscopic quantum self-trapping becomes obvious. The population imbalance exhibits no dynamics within the 
experimental errors and reveals the expected nonzero average (z) ^ 0. Clearly the phase is unbound and winds up over 
time. The error bars in the phase measurements denote statistical errors arising from the uncertainty of the initial population 
imbalance. 



pendulum reaches the top position and continues to ro- 
tate with a non vanishing angular momentum. 

In order to fully characterize the evolution of the sys- 
tem we measure not only the population imbalance but 
also the relative phase of the macroscopic wave functions. 
This is achieved by releasing the Bose-Einstein conden- 
sates from the double- well potential after different evolu- 
tion times. After a time of flight of 5ms in the Josephson 
and 8ms in the self-trapping regime the wave packets in- 
terfere unveiling the relative phase in a direct way since 
the resulting atomic fringes are similar to a double slit 
diffraction pattern. 

In Fig. |21 we present the quantitative analysis of our 
experimental results. The measured fractional popula- 
tion imbalance and the relative phase in the regime of 
Josephson oscillations (z(0) = 0.28(6) < Zc) are shown 
in Fig. Ufa). As expected for a symmetric double- well 
potential the relative population oscillates around its 
mean value (z) = 0. The relative phase of the two 
Bose-Einstein condensates oscillates with a finite am- 
plitude of <f> = 0.5(2)tt around (<f>) = 0. The self- 



trapping regime can be reached by simply increasing the 
initial asymmetry of the double-well potential as indi- 
cated in the schematic diagram in Fig. |2Ib) realizing 
z(O) = 0.62(6) > zc- In this case theory predicts that 
z exhibits only small amplitude oscillations which never 
cross z — i.e. (z) ^ 0. Additionally the relative phase 
4> is unbound and is supposed to wind up in time. In 
Fig. [2Kb) these characteristics of macroscopic quantum 
self-trapping are evident. The population difference does 
not change over time within the experimental errors and 
the phase increases monotonically. The initial deviation 
from the linear time dependence of the phase is due to 
the finite response time of the piezo actuated mirror. 

The experimentally obtained results can be understood 
quantitatively by going beyond the two mode model 
which assumes stationary wave functions in the individ- 
ual wells which is only justified for Ni + N r <C 1000 
atoms ||. Therefore we numeric ally integrate the non- 
polynomial Schrodinger equation [ijj using the indepen- 
dently measured trap parameters and atom numbers. 
The calculations also include the fact that the piezo ac- 
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FIG. 3: Quantum phase-plane portrait for the bosonic 
Josephson junction. In the regime of Josephson oscillations 
the experimental data are represented with filled circles and 
in the self-trapping regime with open circles. The shaded 
region, which indicates the Josephson regime, and the solid 
lines are obtained by solving the coupled differential equ. 
with our specific experimental parameters. The two mode 
model explains the observed z(4>) dependence reasonably in 
both dynamical regimes. The error bars represent the statis- 
tical error and mainly result from the high sensitivity of the 
relative phase on the initial population imbalance especially 
for long evolution times. 

tuated mirror initiating the Josephson dynamics reaches 
its final position only after 7ms. It is remarkable that all 
experimental findings are in excellent quantitative agree- 
ment with our numerical simulation without free param- 
eters. 

The distinction between the two dynamical regimes - 
Josephson tunneling and macroscopic self-trapping - be- 
comes very apparent in the phase-plane portrait of the 
dynamical variables z and <f>. For our experimental situa- 
tion this is shown in Fig. [3] where we compare our results 
with the prediction of the simple two mode model. From 
our experimental observations the critical population im- 
balance can be estimated to zg = 0.50(5). In the frame- 
work of the two mode model [lfj this yields A = 15(3). 
The corresponding solutions of equ. ^ are depicted with 
solid lines. Clearly the basic features of the dynamics are 
well captured by this approach. In the nonlinear Joseph- 
son tunneling regime (z < Zc) the dynamical variables 
follow a closed phase plane trajectory as predicted by 
the simple model. Very recently a variable tunneling 
two mode model has been discussed by Ananikian and 
Bergeman 20] which is quantitative agreement with our 
experimental observations. 

The successful experimental realization of weakly cou- 
pled Bose-Einstein condensates adds a new tool to 
quantum optics with interacting matter waves. It 
opens up new avenues ranging from the generation of 
squeezed atomic states [2]] and entangled number states 



(Schrodinger cat states) |22j to applications such as atom 
interferometry |2^. Moreover the detailed investigation 
of the self-trapping phenomenon could provide a test of 
the validity of the mean field description in atomic gases 
in the strong nonlinear regime |24| . 
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